

NASA CR- 132358 


<!5ASA-CB-1 32358) ANALYSIS OF THE FLOW N74- 14745 

BOOT DELTA BINGS WITH LEADING EDGE 
S5PAEATI0N AT SUPERSONIC SPEEDS (Calspan 

Corp., Buffalo. N.Y.) 43 d HC $4.25 Uoclas 

CSCL 0 1A G3/01 27974 


ANALYSIS OF THH FLOW ABOUT DELTA WINGS WITH 
LEADING EDGE SEPARATION AT SUPERSONIC SPEEDS 


By Joseph P. Ncnni and Chee Tung 


Prepared Under 
Contract No. NAS1-11577 
Calspan Corporation 
Buffalo, New York 14221 



for 


«r 


I 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


ANALYSIS OF THE FLOW ABOUT DELTA WINGS WITH 


LEADING EDGE SEPARATION AT SUPERSONIC SPEEDS 

By Joseph P. Nenni and Chee Tung 
Calspan Corporation 


I. SUMMARY 

A research program was conducted to develop an improved theoretical flow 
model for the flow about sharp edge delta wings with leading-edge separation 
at supersonic speeds. The flow model incorporates a representation of the 
secondary separation region which occurs just inboard of the leading edge on 
such wings and is based on a slender-wing theory whereby the full three- 
dimensional problem is reduced to a quasi two-dimensional problem in the 
cross-flow plane. 

The secondary separation region was modeled by a surface distribution 
of singularities or a linearized type of cavity representation. The primary 
vortex and separation were modeled by a concentrated vortex and cut in the 
cross-flow potential which represents its feeding sheet (in the sense of 
Brown and Michael). The formulation is made determinate by requiring that 
the stream lines in the cavity region be conical rays, by imposing a cavity 
closure condition, by applying a Kutta condition at the wing leading edge, 
and by requiring that the primary vortex and its feeding sheet be force free. 
The flow model reduces to that of Brown and Michael in the limit of zero 
cavity width. 

The cross-flow solutions for the cavity model were obtained, but these 
solutions have physical significance only in a very restricted range of 
angle of attack. The reasons for the failure of the flow model are discussed. 
The analysis is presented so that other interested researchers may critically 
review the work. 

The second order corrections to slender-wing theory for this flow model, 
for nonslendemess and compressibility, have been formally derived using 
matched asymptotic expansions. In the limit of zero cavity width, these re- 
sults reduce to the second-order-theoiy previously obtained by the authors 
using the Brown and Michael cross-flow model. In this case certain improve- 
ments have been implemented in the method of obtaining the surface pressure 
distributions from the velocity potential. The agreement between theory and 
experiment is thereby improved over a limited Mach number range. 


II. INTRODUCTION 

The flow field about low aspect ratio or slender wings that exhibit 
leading-edge separation has been of concern to designers of high-speed 
aircraft and lifting reentry vehicles for some time. The unique features 



of such flow fields are the spiral sheets of vorticity that emanate from 
highly swept leading edges and dominate the flow field on the lee side of 
the wing or vehicle. This separation produces a nonlinear aerodynamic 
behavior of the wing as angle of attack is varied. 

Although the main features of these flow fields have been recognized 
foi almost twenty years, no completely satisfactory theoretical analysis 
has been developed. Polhamus (Reference 1) has developed an accurate ana- 
lytical method utilizing an intuitive leading-edge suction analogy for 
estimating overall forces on slender wings. However methods to predict 
detailed flow quantities such as surface pressures and loadings are needed. 
Moreover, a further understanding of such vortex flows would aid in the 
design of optimum wings for supersonic flight and open up the possibilities 
for control of the flow over such wings. 

The problem is fundamentally non-linear through the boundary conditions 
because the location of the spiral vortex sheets is unknown a priori. However, 
for compressible flow the exact equations of motion are also non-linear. 
Existing analyses for the detailed flow quantities have used slender-wing 
theory whereby the problem is reduced to a quasi two-dimensional incompressible 
problem in the cross-flow plane. The nonlinearity in the boundary conditions 
may then be handled by the use of complex variable theory in the cross-flow 
plane. The most notable efforts in this vein were made by Brown and Michael 
(Reference 2 ), Mangier and Smith (Reference 3), and Smith (Reference 4). These 
analyses differ only in the representation of the spiral sheet in the cross- 
flow plane. Brown and Michael use a concentrated vortex connected to its 
corresponding leading edge by a cut in the cross flow velocity potential. 
Mangier and Smith used a more realistic spiral representation of the vortex 
sheet terminated by a cut and concentrated vortex in the core area of the 
spiral sheets. Smith refined this representation still further by essentially 
considering more wraps of spiral sheet. A refined representation of the core 
area has been considered by Mangier and Weber (Reference 5) . 

Comparison between theory and experiment shows that the existing slendei- 
wing theories greatly overestimate the forces on low apsect ratio wings. (For 
example, the slender wing result is high by almost a factor of two in the 
moderate angle of attack range.) The possible reasons for this discrepancy 
have been attributed to an inadequate cross-flow model and the neglect of 
higher order aspect ratio and Mach number terms in the equations of motion. 

A second order theory was developed by the present authors, in Reference 
6, to account for non-slendemess and compressibility effects. The Brown and 
Michael cross-flow model was used for ease of analysis since «.he gross pre- 
dictions of this model were similar to those of References 3 and 4. The 
consideration of these corrections improved the correlation between 
theory and experiment, but did not remove all discrepancies. It was, however, 
clear from the results of Reference 6 that, in terms of total lift on the 
wing, compressibility and nonslendemess effects were more important than 
refinements in the representation of the primary separation and vortex sheet. 
The essential fault with either flow models of Reference 2 and 3 is that the 



predicted vortex positions are close to the wing surface, which results in 
overprediction of the upper surface suction pressures. It may he concluded 
that a fundamental change in the cross-flow model is required to improve the 
prediction of vortex location and strength. In short, it was felt that neither 
cross-flow model incorporated all of the features that may be physically 
significant to the flow field. 


The major features of the cross-flow as observed from flow-field and 
surface-flow studies (see for example. References 7, 8, and 9) are shown 
sketched in Figure 1. The details of the secondary separation region have 
not been thoroughly explored. In some instances additional flow separations 
and reattachments may be observed close to the upper surface leading edge; 
however, further categorization of this flow region does not appear to be 
warranted at this time. All of these major features except the secondary 
separation region have been considered to some degree of approximation in the 
previous analysis. Hence, it was suspicioned that the secondary separation 
region must be modeled to resolve the discrepancies between theory and experi- 
ment. Moreover, although the secondary separation occurs only over a minor 
portion of the wing span, it can have important effects upon the manner in 
which the flow leaves the wing surface and forms the primary spiral vortex 
sheets . 


A rational, but unsuccessful, attempt to include secondary separation 
in the cross flow model is described in this report. The secondary separation 
was modeled by a surface distribution of singularities and the primary sepa- 
ration in the Brown and Michael sense. In effect the separation is modeled 
as a cavity and in the limit of zero cavity width the flow model reduces to 
that of Reference 2. Unfortunately the mathematical solutions for this flow 
model have physical significance only in a limited angle of attack range near 
zero angle of attack. The reasons for the failure of this model are discussed 
and the results are presented to allow other interested researchers to review 
the analysis critically. 

The cavity flow model was formally incorporated into the previously 
developed second-order theory and in the limit of zero cavity wiHth this 
theory reduces to the results of Reference 6. 

In Reference 6 the surface pressure distributions were obtained from 
the velocity potential by using an approximate form of the isentropic 
pressure-velocity relationship obtained by series expansion. In unpublished 
calculations this expression was found to give spurious results at high 
Mach numbers. The series expansion employed was not valid for high Mach 
numbers. This shortcoming of the previous theory has been partially overcome 
by using the exact isentropic pressure-velocity relationship to calculate 
surface pressures from the velocity potential. This can only be accomplished 
for limited Mach number and angle of attack ranges because, at sufficiently 
high Mach numbers and/or angles of attack, the theory predicts surface 
velocities which exceed vacuum conditions. (This is also a consequence of 
the vorticies being too close to the wing.) In this limited range, however, 
the correlation between theory and experiment is improved. This work is 
described in Section IV of the report. 


3 



III. LIST OF SYMBOLS 


O- 

fk 

c 

C L 

C P 

C pv » C 

■fin) 

G,(S) 

h 

H 

"3 

H + 

l.P. 

* 

PlS) 

Q(s) 

r 

Re 

R.P. 

s,n>* 


ratio of wing span to root chord 

aspect ratio 

cavity width in % plane 

lift coefficient 

normal force coefficient 

pressure coefficient 

pressure coefficient at vacuum conditions 
source strength in cavity region 
portion of solution for $ * 

normalized half span of wing (it is a function of % and i 
times written as fur) to emphasize this dependency) 

function defined in Equation 8 

function defined by Equation 22 

function defined by Equation 24 

a FT 

imaginary part of complex function 
free- stream Mach number 

real coordinate of vortex core in physical plane 

line doublet strength for outer solution 

imaginary coordinate of vortex core in physical plane 

line source strength for outer solution 

distance from 5 axis in S.qX space ( r = */rf~+ 2 ' ) 

Reynolds number based on root chord 

real part of imaginary function 

strained outer variables 


f*- 


s some- 
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S. Y,Z - strained inner variables 
T - normalized cavity thickness 
(J,, - free stream velocity 

W n - complex potential for first inner solution 
W f2t - complex potential for second inner solution 

- Cartesian coordinates normalized by one-half 
the root chord 

a - angle of attic.. 

p - 7 

y - normalized vortex strength (r/h) 
r - vortex strength 
8 - parameter defined in Equation 7 

y - ratio of specific heats 
e - angle of attack parameter ( sin a / ta*> Ji) 

A - real coordinate of vortex in XP lai *® 

A c - inboard edge of cavity in physical plane 

- straining of % coordinate 
cr - complex inner variable, Z + 

(T, - position of vortex in cr plane 

A - delta wing semi-apex angle 

- imaginary part of X 


r 

4 

i* ** 
x n* x »n ,r ui 

X 


X, 




imaginary coordinate of vortex in \ plane 
velocity potential 

inner solution for velocity potential 
complex transformation of <r plane 
position of vortex in x plane 
decomposition of * 


S 






ttl* I 










IV. OUTLINE OF MATHEMATICAL FORMULATION 

The general perturbation method presented in Reference 6 for determining 
the velocity potential is employed again and will only be outlined briefly. 

The coordinate system used is shown in Figure 2. This system is a body 
axis system and the wing leading edge is given by fy = a.h(x). The physical 
spatial variables are normalized by one-half the root chord such that the wing 
apex is located at x ~-f and the trailing edge at x = +f . Then, for a delta 
wing, h--\(1 + x) . 

For infinitely thin wings the velocity potential is expanded as 

§ — | x cos oc +■ sen oc. f sen 1 oc 

and the independent spatial variables may be strained in the typical P.L.K. 
fashion as (see Reference 6) 


y s 7 

x -&*-£, sin at 

The straining in the axial direction, , is determined by the principle 
that "the higher order solutions shall be no more singular than the first 
order solution." This expansion procedure decomposes the original nonlinear 
problem into a sequence of linear poblems. The governing partial differential 
equations and boundary conditions on the are given by Equations 7, 10, 

12, and 15 in Reference 6. 

For wings of small aspect ratio the are determined by i.. a method of 
matched asymptotic expansions. The small parameter used is a,, the ratio of 
wing span to root chord. The outer expansion for valid far from the wing 
and for small a, is of the form 


iff ~ 7 * cl * cl 4> f + o(a, ¥ ) ( 1 ) 

~ 4/ D 

where^ ff, = land ff t j is the potential for a line source distribution 

and 4^, is the potential f or a line doublet distribution. The doublet distri- 
bution is required for lifting flows and the source distribution is required 
because the cavity representation for the secondary separation region 
effectively gives the wing thickness. The appropriate forms are 



J_ P S ' Ar Q(s 9 )ds m 

2rr *4 [(S "S 0 )*-/5‘*r 


( 2 ) 


o 1 n r s ' flr p(s<>){s-s e )ds 0 

m «r r /, [<*-*) ] ,,M 
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where 


"fX* 


J ■ 






% ■ 


m 4 l 

r * r? + $ 


fi* ' - / 


The source strength Q(S 0 ) and the doublet strength P(S a ) will be determined 
later by the matching conditions between inner and outer solutions. The 
inner problem is obtained by stretching the coordinates normal to the wing 
to regain the detail lost near the wing. The inner variables are 

Y - % /<*■ 

l C/a. 

S -- S 

Then the inner solution valid near the wing for small values of a. is of the 
form 


**/ - a if, ( y. Z ; S ) * a* in, ( Y. i ; s) + d/in,* ( Y,l,s ) 

v- o(aX t a ¥ Jm.a.) 


(4) 


The partial differential equations for the are given by Equations 26, 27, 
and 28 of Reference 6. The important point to note is that these equations 
are two dimensional in Y and ? with 5 appearing only as a parameter. The 
first inner solution satisfies LaPlace's equation in the plane transverse 
to the wing and the appropriate form of surface flow tangency. This first 
inner solution is the standard "cross-flow problem" of slender-wing theory 
and may be solved by conformal mapping techniques of complex variable theory. 
The first inner solution may be expressed as 

#/, - K.P. yVff(o-) + G, (S) 


where <r is the complex variable £ and \V„ is referred to as the complex 
velocity potential. The G,(S) term is required because of the source-like 
behavior of the cavity, and it is determined by matching with an outer 
solution. 

A brief review of the experimental observation of the cross-flow and the 
detailed formulation of the cross-flow model used is presented in the next 
section. 


V. THE CROSS-FLOW MODEL 

Review of Experimental Observations of 
the Cross Flow 

The tentative concept of the actual cross-flow model for slender wings is 
shown in Figure 1, deduced largely from the vapor screen studies of References 
7 and 10, the smoke flow studies of References 8 and 11, and the surface flow 
studies of References 13 and 14. (These references have been selected because 
of the quality of the photographic material contained in them rather than on 
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chronological order or originality.) These flow features are present 
regardless of the Mach number as long as the leading edge is subsonic. 

There are some systematic differences evident in the cross-flow model at 
subsonic and supersonic speeds. The vapor screen studies of Reference 8 
indicate that at supersonic speeds the vortex pattern is spread over a 
greater spanwise extent of the wing and lies closer to the wing surface than 
at subsonic speeds. Comparison of surface pressure distributions (as between 
References 10 and 12) also confirms the observation that the vortex pattern 
is more diffused at supersonic speeds. 

Details of the secondary separation region are dependent upon Reynolds 
number. The spanwise flow evidently separates because it cannot negotiate 
the pressure rise between the negative pressure peak (underneath the vortex 
core) and the leading edge. Available experimental data indicate that the 
primary Reynolds number effect is through boundary layer transition and its 
proximity to the secondary separation line. There is some experimental 
evidence which indicates that for the higher aspect-ratios, reattachment 
occurs between the secondary separation region and the leading edge. 

An example of Reynolds number effects on the surface pressure distribu- 
tions may be seen in Figure 3 which shows some subsonic data from Reference 11. 
As Reynolds number increases, the magnitude of the upper surface pressure 
peak increases, and the lateral extent of the vortex pattern decreases. 
Apparently, after the boundary layer at the secondary separation point be- 
comes fully turbulent, the Reynolds number dependence becomes very weak. The 
spanwise pressure distributions also indicate that at the lower Reynolds 
numbers the secondary separation region contains several discrete longitudinal 
vortices (as indicated by the variations in pressure) while at higher Reynolds 
numbers the secondary separation is more like a cavity (as indicated by the 
uniformity of the pressure). Although the pressure distributions are markedly 
affected by Reynolds number, the normal force and pitching moment are 
independent of Reynolds number. (The wing root bending moments which have 
not been measured, however, must be affected.) 

It must, however, be concluded that the existence of any secondary 
separation region is a significant factor in the observe-d differences between 
theoretical predicted forces and experimentally determined forces. Since the 
theoretical values are higher, the foregoing discussion would indicate that 
the secondary separation region of slender wings is an area where boundary 
layer control devices may be used to increase the lift on such wings. Complete 
elimination of the secondary separation region should move the experimental 
results closer to the theoretical results. 

The experimental investigations of Reference 14 indicate that at transonic 
speeds there are similar effects of Reynolds number or. the pressure distribu- 
tion. However, at transonic speeds the normal force and pitching moments 
are also dependent on Reynolds number. These results cast some doubt as to 
whether the existing supersonic data are free of Reynolds number effects. 

More experimental work is needed in the supersonic speed regime to investigate 
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Reynolds number effects. Moreover, additional detailed flow field and 
surface pressure measurements near the leading edge are needed to investigate 
the character of the secondary separation region. 


Mathematical Flow Model 

The secondary separation region is properly modeled by a volume distri- 
bution of singularities over the separation region. However, since this 
representation is mathematically intractable at present, a simpler represen- 
tation was sought. Since the height of the secondary separation is small 
compared to the other wing or flow field dimensions, it appeared permissible 
to collapse the volume singularity distribution to a surface singularity 
distribution. Indeed, inspection of the vapor screen studies of Reference 7 
would suggest that at low angles of attack even the primary separation might 
also be modeled by a surface singularity distribution. (This type of flow 
model was considered further and will be discussed later in Appendix A.) Such 
a flow model would then be generically related to that used in Reference IS 
to analyze the flow over a partially cavitated, low-aspect-ratio hydrofoil. 
Hence, the terminology of cavity has been used for this surface singularity 
distribution. For reasons discussed pre iously the Brown and Michael represen- 
tation of the primary separation was considered adequate for present purposes. 

Inclusion of the secondary separation in the flow model requires the 
development of additional boundary conditions to make the problem determinate. 
Since the condition of flow tangency at the wing surface must be abandoned in 
the cavity region, one of the other flow quantities must be prescribed on 
the cavity. In the hydrofoil case, Reference 15, the pressure is specified 
as constant along the cavity, but tnis is not appropriate for aerodynamic 
problems. The required constraint for the present flow model was deduced 
from inspection of surface flow patterns, primarily those shown in Reference 
9. Inspection of these photos shows that in the secondary separation region 
the surface streamlines are conical rays emanating from the wing apex and 
this was the condition used, being equivalent to making the transverse velocity 
parallel to the wing in the cross flow plane equal to u m f/s. Since the wing 
surface in the cross-flow plane also has a transverse velocity of £C/S, this 
condition implies that in the cavity region of the cross-flow plane there 
is no relative transverse velocity between the flow and the wing. Then, 
according to the definitions of separation in three-dimensional flow given 
by Maskell in Reference 16, the inboard edge of the cavity is automatically 
a separation line. In total, then, the conical streamline requirement in 
the cavity region appears to be consistent with all of the physical require- 
ments for modeling the secondary separation region. 

One additional constraint is required to make the problem determinate. 

The condition selected was that of cavity closure. To the inviscid flow 
outside the cavity, the cavity appears as a local thickening of the wing. The 
cavity closure condition is then the requirement that the cavity has no 
thickness at its edges, implying that the flow that leaves the wing surface 
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at the inboard edge of the cavity reattaches to the wing at the leading edge 
before departing from the wing surface and becoming entrained in the primary 
separation flow pattern. 

The resulting cross-flow model and applicable bounds' • value 'n the 
velocity potential are shown in Figure 4. The subscript notation has ht>>' 
used to indicate partial differentiation with respect to_thc appropriate vari- 
able. The concentrated vortices are located at <r\ and-?; ind are connected 
to their respective leading edges by a cut in the cross-flow potential. The 
cut then represents a planar vortex sheet in the three-dimensional flow with 
the axis of the vorticity vector in the cross-flow plane. In this fashion the 
cuts serve as the feeding sheets for the concentrated vortices. The cavity 
extends between \ and h . 

This boundary value problem is most easily handled by using conformal 
transformations. The transformation chosen was 


K * (5) 

The branches of the square root are chosen such that the upper wing surface is 
mapped into a line segment along the positive imaginary axis and the lower 
wing surface is mapped into a line segment along the negative imaginary axis. 
The flow in the * plane is shown in Figure 5. Any analytic function with the 
proper symmetry in this plane automatically satisfies the condition of no 
flow through the wing. The cavity is represented by a source distribution 
in this plane that extends between the origin and ic. The vortices are 
mapped to X, and-*,. The appropriate complex potential is then given by 


Wnti) - -i< f(q 0 ) Jn(x-Lr] 0 )4ri 0 - X, (<0 


Whereas the first term gives the flow about the wing without any additional 
singularities in the flow field, the second term gives the flow due to a 
distribution of sources of local strength -f(n t ) along the cavity, and the 
last term gives the flow due to the concentrated vortices . The source 
strength is given by the requirement for conical streamlines in the cavity 
area. In terms of the velocity potential this is 


K - « 


on Y - 0 , 




(7) 


{ - -A 


Utilizing the relation that f. 1 . * R.P. and Equation 6, the following 




d.ar 


integral equation results for/(») . 


(«) 
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The inversion of this equation can be obtained from Reference 17 or 18 as 


+(n) 



i , 

%-c *7 ^ 


L j 

f c - 7 

\ „ 

. r J 

k- 7 

. 1 r i 

IJL 

* 1 

7 

i 

7 1 


1 K^t-C 


The solution to Equation 8 is not unique in that a term of the form co»st /-/r[Uj -c) 
may lie added to f(n) and still satisfy the equation. This term, however, would 
introduce undesired singular behavior at the edges of the cavity and therefore 
has been discarded. 

The leading edge Kutta condition or requirement for smooth outflow at the 
wing edges is 

iiL'l *0 

dx. I V , (10) 


which implies that 


f(o) = 0 


or 



r 

[ / 

nr + j_ 

1 Xf 

ITT 

1 */ 1 

f X r iC X, ) 

n \,+ic 


0 


( 11 ) 


In the limit c -* o this condition reduces to the Kutta condition for the Brown 
and Michael flow model. 

For conical flow the outer boundary of the cavity should be expressible 
as 

r* 


Then the requirement that the outer edge of the cavity is a stream surface is 
that 


(<>•*) - (* T )s * 0 
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Using the conical flow assumptions this may be expressed as 


hr = if' ft, (0 t £ o ) * rrzj ? f(q 0 ) 

* *o ° 'It 


Then the requirement for cavity c' »ure is T * 0 at c - O 
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- * 

'O h -q. (o 


or 




ir_ 

’ Zrr 


Ki f i IEHE rl L I fEU 
^-icVh 2 + t*L» K, J 2h(h + iK,)\j h 


2h(h-i\ ] ) 




( 12 ) 
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’JE!.,' 

f 



IV* C Uh* 

1 *, 

2b.(h + L \, ) 

1 h 

2h(h-ii t ) 

ii~ 1 



- O 


In the limit C-*0this equation becomes an identity. 

The force balance on the concentrated vortex and feeding sheet (cut in 
complex potential) is given by 


e d*;, 

2 do- 


d crj 

ds 


dr 

ds 




(13) 


dw* dw„ 

where --j— is evaluated at a~ t less the singular contribution from the vortex 
locatedTat o~ f .For the conical flow F and o~, scale with h and Equation 13 becomes 




r 1 ire t _ 

1-1 ♦ - ( 2<r , 

Inh 


i) - o 


(14) 
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where 



* rc ' £_ r*i | l \,y r _ I 

tirfX' + X') lx, [7 l, 1 %, , 


In the limit c —0 this equation reduces to the corresponding force balance for 
the Brown and Michael flow model. 

Equations 11, 12, and 14 then provide sufficient conditions to determine 
the unknown parameters c , r , and <77 . This system of equations must be solved 
numerically. The details of how this was accomplished are given in Appendix B, 


Results and Conclusions 

The system of equations for the unknown parameter C, r and <7- is non- 
linear and hence the number of roots is unknown. It may be seen by inspection 
of the equations that Equations 11 and 14 reduce to the corresponding syst< m 
of equation for the Brown and Michael cross-flow model in the limit c-»0. 

Also, Equation 12 becomes an identity in this limit. Therefore, the Brown 
and Michael values for the vortex parameters and C - 0 will satisfy the system 
of equations. This is not the desired solution and will not be discussed 
further. 

The only solutions with nonzero c that have been found are shown in 
Figures 6, 7, 8, and 9. These roots have physical significance only below 
€ • 0.206 and hence the figures are restricted to this range. As seen in 
Figure 6, above Z = 0.206, c becomes negative which implies that the cavity 
shifts to the wing lower surface. This clearly contradicts experimental 
observations. Furthermore, in the range of e where c is positive, c is 
very small and the corresponsing vortex parameters are within one half of 
one percent of the corresponding Brown and Michael values. (Therefore these 
values have not been included in the figures.) Wing surface pressure distri- 
butions for? = 0.1 (near the maximum cavity width) are shown in Figure 10 
for the Brown and Michael cross -flow model and Figure 11 for the cavity 
cross-flow model. Although the cavity does not significantly change the 
vortex position or strength, the upper surface suction peak is slightly 
reduced, producing a slightly lower normal force coefficient. Although these 
changes are in the right direction, the magnitude is inconsequential and 
therefore this specific cavity cross-flow model has little application to the 
slender wing problem. 
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The failure of this flow model must he due to either the conical 
streamline condition imposed upon the cavil) or the cavity closure condition. 

In the present flow model it was assumed that the cavity closed at the wing 
leading edge. It is possible that the cavity may close inboard of the 
leading edge (some argument may be made for this assumption by inspection 
of surface flow studies of Reference 13) or may not close at all. The present 
experimental observations are insufficient to answer this question con- 
clusively. 

A flow model in which it is assumed that the cavity closes inboard of 
the leading edge contains an additional unknown parameter, namely the location 
of the outboard edge of the cavity. It appears that such a model could be 
made determinate by eliminating the singular behavior at the outboard cavity 
edge, but this flow model has not been investigated in detail. 

An alternative condition to replace the conical streamline condition on 
the cavity might be deduced from that used in Reference 15. The cavity 
cross-flow model used in Reference 15 for the cavitating low-aspect-ratio 
hydrofoil problem evidently produces realistic results for that problem. In 
that analysis a constant pressure boundary condition is imposed on the cavity. 
The specific form used was that the spanwise velocity component on the cavity 
was a constant, which is tantamount to assuming that the vorticity in the 
cavity is constant. This interpretation of the boundary condition used might 
be plausible for an aerodynamic problem. It is difficult to assess before- 
hand how this would affect the solution for the model presently used. However, 
the flow model used in Reference 15 results in singularities at both edges of 
the cavity. The leading edge singularity may be removed by adding the primary 
vortex to the flow model, but it is not immediately evident how the singularity 
at the inboard edge of the cavity may be removed. 

There are then several plausible modifications to the present flow model 
that might be considered but essentially on a trial and error basis. However, 
there appears to be a more systematic way to address the problem. Thus far 
the relevant theories have been using potential flow in an attempt to model 
a flow in which viscosity plays an important role. (The primary and secondary 
separation are directly attributable to the presence of viscosity.) Any 
potential, flow model which attempts to model these features is inherently 
indeterminate. The attempts to resolve these indeterminacies have been to 
deduce the proper boundary conditions by studying experimental observation 
in which viscosity is inherent. Thus far, however, the proper boundary 
conditions have not been completely identified. Moreover it appears that 
the surest way to significant progress on sol 1 - ing these problems would be to 
consider a flow model in which the viscous terms are present in the equations 
of motion thus eliminating the need for adu.'tional boundary conditions. 


VI . SECOND ORDER THEORY FOR CAVITY FLOW MODEL 

The previously described cavity flow model has been formally incorporated 
into the second order theory developed in Reference 6. Recall that 


1 * 
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( 15 ) 


K ~ ? - * % - 


~ A $n * a * *<n - K 


12 % 


( 16 ) 


where 


K = R.P. W„ + &,( s) 


and I4^ r is given by Equation 6, **ij is given by Equation 2, and 4>^ is given 
by Equation 3. Now, using the asymptotic matching principle and equating the 
outer expansion of Equation 16 to the inner expansion of Equation IS results 
in 


and 


Q(s) = ( C f(f\ 0 ) at no ( 1 ?) 

¥r £ {7 ' £ < v <- > -/*> ^ } (i8) 

G,(s) = Cffs) (19) 


£• 

Now $ tZt satisfies the inhomogeneous partial differential equation given by 
Equation 27 of Reference 6 and boundary condition given by Equation 30 
thereof. The method of solution is also the same as that of Reference 6. 
$t*f is composed of a particular part and a homogeneous part and may be 
expressed as 


W 12} (cr) + H 3 $;] 


( 20 ) 
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The unknown portion of the second term of Equation 19 is determined by 
matching of inner and outer solutions to be 
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The JV integral must, in general, be evaluated numerically. The coefficient 
of the logarithmic term in Equation 16, § l satisfies LaPlace's equation and 
may be expressed as 

fL ■ ", K W 

The function is again determined by matching of inner and outer expansions 
to be 


* 


(24) 


To calculate wing surface pressures it is not necessary to determine any 
more of the higher order solutions (see Reference 6). The inner solution for 
the straining is, however, required and this is the same as that determined 
in Reference 6, i.e.. 


£ = -"‘K + r 


where F vanishes on the wing. 


In Reference 6 the surface pressures were calculated from the velocity 
potential by using an approximate series expansion form of the isentropic 
pressure-velocity relationship. This was found to give spurious results at 
high Mach numbers where the series expansion broke down. In the present 
work the exact expression was used. Then on the wing surface the isentropic 
pressure-velocity relationship gives the following expression for the pressure 
coefficient 


C p =■ j[/ ~ ^j---M*^2stnacCOsottff l& + sin 2 a |V,* (1 + 2rt*co$ac ) 

*< f *H)r / - , j 


(25) 


where the quantities in this expression are evaluated on the wing surface and 
are given by the following expressions. 

On the upper surface 


* (f + 2a. 2 Hj +■ 2 cl^oo. 
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and 
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where if * j f-/*and ? is the fraction of half span. 
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where 
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On the lower wing surface 
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where 
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On the lower surface if c + 0 
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These expressions have been used to calculate the surface pressures for 
a case in the e range where CfO for the cavity cross-flow model. These 
results are shown in Figure 11. The limiting forms of the expressions when 
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c going to zero were also used to calculate the corresponding pressure distri- 
bution for the Brown and Michael flow model. These results are shown in 
Figure 10. 

A comparison of the two pressure-velocity relationships is shown in 
Figure 12. It is seen that the exact isentropic relationship is somewhat 
better than the approximate form; however, the upper surface suction peak 
is still overpredicted. The exact expression proves to have a limited range 
of applicability. At sufficiently high Mach numbers and/or angle of attack 
the theory predicts surface velocities that exceed those corresponding to 
vacuum conditions (because the primary vortex is too close to the wing surface). 
A comparison of the lift predicted by the two forms of the pressure-velocity 
relationship provides a better estimate of lift in tlv regions where the 
calculations can be performed (Figure 1’ 


VII. CONCLUSIONS 

A new flow model for the prediction of the flow about delta wings with 
leading-edge separation was investigated. The flow model included a 
representation of the secondary separation region that has been experimentally 
observed on such wings. Secondary separation was modeled as a cavity. The 
resulting solution for the flow model proves to have only a very limited range 
(near zero angle of attack) of physical significance. Thus, these results 
cast suspicion upon the boundary conditions that were used to make the model 
determinate. These boundary conditions consisted of a cavity closure condition 
and a requirement for conical streamlines in the cavity region. These assump- 
tions are plausible, based upon the available experimental evidence, but 
evidently one or both are incorrect. It appears that the correct choice of 
boundary conditions can be developed only by considering a flow mouel that 
retains viscous terms in the equations of motion. 
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Appendix A 

CROSS FLOW MODEL WITHOUT CONCENTRATED VORTICES 


As mentioned in Section V of this report a flow model was investigated 
in which the primary separation was modeled by a surface similarity distri- 
bution. 

The appropriate complex velocity potential is 


w n s 


(A 1 ) 


The conical streamline requirement on the cavity gives 



With this flow model there must be a singularity at either or both edges of 
the cavity depending on the specification of k, . The choice fc = ~ ( 1 + L& ) 
meets a Kutta condition at the leading edge and produces a quarter roo^ 
singularity at A c . The resulting slender-wing normal force is then given by 

It is seen from this expression that the cavity produces a non-linear contri- 
bution to the lift which is opposite in sign of that of experimental 
observation. 

The choice k • 0 produces no singularity at Ac, however results in a 
three-quarter root singularity at the leading edge, (The physical interpre- 
tation of this leading edge singularity is questionable.) The corresponding 
slender-wing normal force coefficient is given by 


'N 


ir 4? f . e 2 / c \ c 2 ) 

— s.ncc I -_j 


(A4) 


It is possible for the nonlinear contribution to be of the proper sign, and 
hence further analysis was conducted. 

In these flow models C is determined by the requirement for cavity closure. 
Mathematically this is r e /(q g ) 

i /T*^r * * 
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or 
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This equation only has solutions for fc>23.3 This large an e value is well out 
of the ranges of practical value. 

Hence it was concluded that cross flow models of this type were of no 
use for the slender wing problem. 
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Appendix B 

NUMERICAL SOLUTION OF THE CROSS-FLOW EQUATIONS 


The system of equations to be solved for the unknown parameters c , T , cr t 
is composed of Equation 11, 12, and 14. The basic method employed was a 
Newton-Raphson iteration scheme. First, the equations are normal ized In- 
dividing all lengths by h and by letting y • r/h . Then y is eleiminated from 
the system by using Equation 11, viz. 


7 


ir( 1- 


Sc \ 

t) 




(Bl) 


and the following functions a-c defined. 


and 


CF12 - (2 t r. -f) 


r 
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* i - 


_L— 

2rr <, + \ z 


(B2) 


FI -- 

R.P. CF tZ 

(B3) 

F2 = 

I P. CFtZ 

(84} 


From Equation 7 




(B5J 


and 


F3 - f‘3/c* 


(B6) 


The values of and c which simultaneously zero the F t - functions are then the 
desired solution to the problem. The Newton procedure is outlined below. 

An initial guess at the solution at a selected 7 value is made, say, \ (t) 
and c(l) , Further let X, * tZ.P. X* * IP-X^ and X tn » c . Then it is assumed 
that if these guesses are close to the solution that 
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F i * F i^ 
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then the increments in the variables to reach the solution arc fount 1 by solving 
the linear system of equations. 


or 


where 


"* - -w 


■' ■ 1 . 2 . 3 fB8J 


- F 1(1) 
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3 F3 

9X, 

3x t 

9x 3 


All the partial derivatives are evaluated at There are similar 

expressions for AX Z and AXg. The next guess at the solution is then 
/■(2) m X±(l) + &X±(f)‘ This process is then repeated until all \F-(n\ < 10 
The partial derivatives of the Fj were computed using the finite difference 
approximation. The initial difference was chosen as . 01*^. The difference 
for the successive steps where chosen as the AX;(rj-f) at each step. 


This scheme was programmed for the Calspan./IBM 370 computing machine. 
The resulting solutions for non-zero c and vortex positions are shown in 
Figures 6, 7, 8, and 9. 


Initially the system was used with F'3 as given by Equation B5. This 
system of equations (with the method of solution used) would only converge to 
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solutions for which C* 0 . Analysis of Equation B5 for small c then showed 
that F'3 vanished as c* for sm«ll c , hence Fi was formed and this system 
resulted in the solutions shown. 

The entire system of equations was analyzed for small c because it was 
intuitively thought that there should be a solution with c si.all, at least 
in the lower range. This analysis showed some properties of the system of 
equations that were interesting in their own right and therefore this analysis 
is outlined below. 

Equations 11, 12, and 14 were expanded for small C . The results are as 
follows. The leading edge Kutta condition. Equation 11, becomes 


where 


F, a» * cf * + k * o(c> ) * 0 
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The closure condition. Equation 12, becomes 

I F 'ari * i c * f < + c \iz F, »* * )+ o(cU) * 0 


(*9) 


(BJO) 


The force balance on the primary vortex and feeding sheet, Equation 14, 
becomes 


c Fj + *■ o ( c*) 9 O 


(Bin 


where 






Ii JjL t . j 1 fL + ill l 

Z 1 Ztr + X, } 2 rr <r t 
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Now it is recognized that fr i e „ - 0 and are respectively the leading 

edge Kutta condition and the force balance on primary vortex and its feeding 
sheet for the Brown and Michael cross-flow model. Hence, as c -*■ 0 Equation Blr 
vanishes as c z . This explains why the original system of equations would onl; 
converge to solutions with c=<7 . 

Furthermore c/Z may be factored out of Equation BIO resulting in 




cf, + c 




I Bid) 


Then Equation B9, Bll and B12 form the approximate equation for the cross- flow 
model for small c . 


Comparison of the approximate closure condition. Equation B12 and the 
approximate Kutta condition. Equation B9 shows that to zeroth order in 
that these conditions are not independent. Because of the special form of 
Equations B9 and B12 they may be conbined to give 


c 



s * z * 7T FlBM 


(B L 3 ) 


so that c may be eliminated from the system of equations. Then Equation B9 
and Bll constitute a sufficient system to determine y and \ f . This approximate 
system must also be solved numerically and this was done in a fashion similar 
to the full equations. The solutions to the approximate system were found to 

have the same general behavior as the solutions to the full system and thus 
reinforced the results previously given for the full system. 


S' 
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Figure 1 OBSERVED CROSS FLOW (SHADED AREAS INDICATE REGIONS OF 
SIGNIFICANT TOTAL HEAD LOSS) 
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Figure 3 EFFECT OF REYNOLDS NUMBER ON UPPER SURFACE PRESSURE 
DISTRIBUTION FROM REF. 13 
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Figure 4 CROSS FLOW MODEL WITH CAVITY IN <T PLANE 









Z, FRACTION OF SEMI-SPAN 


Figure 10 WING SURFACE PRESSURE DISTRIBUTION SECOND-ORDER THEORY WITH 
BROWN AND MICHAEL CROSS-FLOW MODEL, M - 1.5, JR - 1.0, l - 0.1 
( a - 1.43®), C N - 0.0503 
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Z, FRACTION OF SEMI-SPAN 

Figure 11 WING SURFACE PRESSURE DISTRIBUTION SECOND-ORDER THEORY CAVITY 
CROSS-FLOW MODEL, M - 1.6, 4* - 1.0, € « 0.1 ( a - 1.43°). X c * 0.9971, 
C N - 0.0485 
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SURFACE PRESSURE COEFFICIENT. Cp 



LIFT COEFFICIENT, C 




ANGLE OF ATTACK a, DEC 

Figure 13 COMPARISON OF THEORETICAL LIFT COEFFICIENTS USING EXACT 

AND APPROXIMATE ISENTROPIC PRESSURE-VELOCITY RELATIONSHIP 
A * 0.25 
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